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Abstract 

We give a review of the exact renormalization group (ERG) approach and illustrate its applications 
in scalar and fermionic theories. The derivative expansion and approximations based on the derivative 
expansion with further truncation in the number of fields (mixed approximation) are discussed. We 
analyse the mixed approximation for a three-dimensional scalar theory and show that it is less effective 
than the pure derivative expansion. For pure fermionic theories analytical solutions for the pure 
derivative expansion and mixed approximation in the limit N — » oo, where N is the number of 
fermionic species, are found. For finite N a few series of fixed point solutions with their anomalous 
dimensions and critical exponents are computed numerically. We argue that one of the fermionic 
solutions can be identified with that of Dashen and Frishman, whereas the others seem to be new 
ones. The issues of spurious solutions and scheme dependence of the results are discussed. 
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1 Introduction 



One of the methods, which, in principle, is capable of handling non-perturbatively many problems of 
quantum field theory (such as calculation of exact renormalization flows, explicit derivation of low-energy 
effective Lagrangians from high-energy ones, extension of the c-theorem to higher dimensions, etc.) is the 
exact renormalization group (ERG, hereafter). Originally developed by Wilson in his seminal articles in 
the early seventies |jj (see Ref. ^ for a classical review), it has recently attracted much attention. 

In this article we are going to consider a version of the ERG approach based on a functional differential 
equation. It determines how the action changes as momentum modes are integrated out in lowering down 
a characteristic momentum scale (cutoff) A while keeping the S'-matrix elements invariant. This is the 
so-called ERG equation. Using it for the most general action, consistent with the symmetries of the 
model under consideration, the complete set of /3-functions (or How equations) can be derived, and from 
them the location of fixed points, critical exponents, behaviour of the renormalization flows, low-energy 
effective action, etc. can be computed in principle. However, with known techniques it is not possible 
to handle such general actions for physically interesting models. One needs to choose an approximation 
which, on one hand, reduces sufficiently the complexity of the problem and, on the other hand, captures 
the essential features of the physical system under consideration. Approximations most commonly used 
in the literature are the local potential approximation || - , polynomial approximation for the effective 
potential, which can be regarded as further simplification of the local potential approximation, ||], Q 
and the derivative expansion ]l0| - |14| . Discussion of their main features and some results obtained by 
these techniques are given in Sect. 2.3. For calculations we will use the derivative expansion in the case of 
fermions and the derivative expansion combined with expansion in powers of fields (mixed approximation) 
for fermionic and scalar theories. 

In this article we will consider the Polchinski type ERG equation [jl5| with its derivation being dis- 
cussed along the lines of Ref. jL6) . It is worth mentioning that there are other types of ERG equations 
used in quantum field theory calculations. Among them is the one of Ref. ||, based on the sharp cutoff 
realization of the Wilson idea by Wegner and Houghton JItJ , the average effective action equation ]l8| , a 
version of the ERG for the effective action with a sharp cutoff of Ref. |1S|| , to name a few. 

Although firstly used for studies of critical phenomena in condensed matter problems |^ti|] , now the 



scope of the ERG approach covers many other areas 21 , including field theory, as it was demonstrated 
by Polchinski in an elegant paper where he proved the perturbative renormalizability of A</> 4 in a quite 
simple way |jl5|] . Most of the work was devoted to scalar theories, both to their general studies |lE{| , p2f 
and to computation of fixed points and the critical exponents || - || , || - Q , [jll) - Jn| . Much work 
has been done on gauge theories, see Refs. [^3| for examples. However, a satisfactory formulation of 
the approach is lacking in this case, the main difficulty being to maintain the gauge invariance along the 
renormalization flow of the regularized action. It is worth mentioning studies of phase transitions at finite 
temperature |^4|, of bound states |2q] , p6| and in relativistic cosmology 27 within the ERG approach. 
Investigations of theories with fermions, which have been done so far, are very restricted. In p8| the 
ERG approach was used for proofs of the renormalizability and convergence of the perturbative expansion 
of the two-dimensional Gross-Neveu model. Discussions of physical effects in many of the articles are 
limited to perturbative expansions p9[ or models with fermions interacting with scalars either through 
Yukawa coupling p0[ or through coupling of the fermion bilinears to a scalar potential and then studying 
the fixed points and renormalization flows for this potential [Q , |2(| , |3lJ . The ERG formalism for pure 
fermionic theories, in which fermionic degrees of freedom are treated directly and non-perturbatively, has 
been developed recently in Ref. |3^ ]. Fixed point solutions and the critical exponents were also analysed 
there. 

The purpose of this article is to give an account of basic ideas of the approach and main steps of the 
derivation of the ERG equation (for the scalar and fermionic theories) and to describe main techniques 
used for non-perturbative calculations of fixed points and critical exponents. We will also give a review of 
main results of such calculations and illustrate the application of the techniques in two concrete models. 
One of them is a three-dimensional model of one-component scalar field with ^-symmetry. It possesses 
the non-trivial fixed point (Wilson fixed point), and was analysed by many authors (see j|, Q - 
H, JllJ - JfJ]). The other example is a two-dimensional Gross-Neveu type theory which seems to be a 
perfect site for such studies: its phase diagram is non-trivial and, on the other hand, the two-dimensional 
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Clifford algebra is the simplest one and generated by the Pauli matrices. In the analysis of this model we 
will follow mainly Ref. p2[ , but also some new results on the pure derivative expansion will be presented. 

We would like to mention that pure fermionic theories, analysed in this article, are widely used in 
description of fundamental interactions. There are some interesting phenomenological models of this 



kind, like the celebrated Fermi theory of weak interactions 1 33 , models of resonance physics [B4j based 



on extensions of the Nambu-Jona-Lasinio action |35[, models explaining the symmetry-breaking sector 
of the electroweak theory, especially in connection with technicolor theories J3(|, etc. Also we would like 
to remind that fermions are not always easily manageable by non-perturbative methods, e.g. in lattice 
formulation. On the contrary, as we will show, there are no principal difficulties in treating fermions 
within the ERG approach. Moreover, once truncated, the ERG equations for fermions and bosons look 
very similarly, thus making possible applications of techniques, already known in bosonic theories, in the 
case of fermions. 

Moreover, the two-dimensional Gross-Neveu model E^J , which contains N species of fermions trans- 
forming under a global representation of the unitary group U(N), is interesting in its own. It is asymptot- 
ically free and renormalizable within perturbation theory and also within the 1/iV expansion. However, 
none of these approximations is capable to find any non-trivial fixed point for d — 2. An interesting 
modification leads to the so-called chiral Gross-Neveu model |3lj , which is chosen to have the additional 
symmetry of the U(l) chiral group. As in the previous case no fixed points, besides the Gaussian one, 
can be found within the l/N approximation. However, using non-perturbative methods based on current 
algebra and conformal techniques Dashen and Frishman in Ref. ]38| | found two critical curves in the 
space of couplings for which the theory is scale invariant. One of the lines corresponds to the abelian 
Thirring model, whereas the other one is truly non-trivial and does not pass through the origin. A very 
remarkable fact is that this result is exact and is not given by a zero of a /3-function, neither in the 
perturbation theory nor in the large N expansion. For this two continuum sets of critical theories the 
value of the coupling constant associated to the SU (N) degrees of freedom is fixed to be equal to zero or 
47r/(AT+l) respectively, while the coupling associated to the abelian degrees of freedom is arbitrary. More 
recently, using bosonization, current algebra and conformal techniques, other non-trivial fixed points in 



two-dimensional fermionic models were found |39 



The article is organized as follows. In Sect. 2 we derive the ERG equation for scalar theories and 
study the mixed approximation for a three-dimensional model. Results of computation of the fixed 
points and critical exponents are compared with those obtained by the pure derivative expansion and the 
polynomial approximation in earlier works. In Sect. 3 we derive the ERG equation for pure fermionic 
theories for arbitrary number of spacetime dimensions. Sect. 4 is devoted to the construction of the action. 
The chiral Gross-Neveu type model, we are going to study, is defined through its symmetries, and the 
approximations, we are going to use, are explained. In Sect. 5 we work out the pure derivative expansion 
and obtain analytically two fixed point solutions. The calculation of the /3-functions, fixed points and 
critical exponents within the mixed approximation in the fermionic model is covered by Sect. 6. Sect. 7 
contains discussion of the results and some conclusions. 



2 Scalar theory 

In this section we explain the main points of the derivation of the ERG equation for the scalar case 
following the lines of Ref. [j^. We also present an example of calculation of the fixed point and the 
critical exponents in a scalar ^-symmetric model in three (Euclidean) dimensions in order to illustrate 
how the whole machinery works. 

As usual in this type of problems, there are three important steps. The first one is to derive an ERG 
equation which governs the behavior of the action as we integrate out modes, that is, how the action of 
our effective theory has to be modified when we vary the characteristic scale (cutoff) A. This is achieved 
by requiring the independence of the S-matrix elements on A. However, it is more convenient to impose 
a stronger condition, namely the independence of Green's functions on A. The second step is to choose 
the appropriate "space of local interactions" . This amounts to define properly the theory and a sensible 
approximation of the general action. The last step is to compute the /3-functions, which characterize the 
RG flow of the theory, and to obtain physically relevant information from it. We may be interested in 
global properties of the flow, like the number of the fixed points it contains, or just local ones, like the 
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behavior of the flow in the vicinity of them. The latter is characterized by the critical exponents, which 
are, moreover, universal quantities. 



2.1 ERG equation 

To derive an ERG equation we impose the condition that the generating functional of Green's functions 

Z[J] = J Z^exp (-S [</>; A] + J JpQ^-p + /aJ 

is independent of some cutoff A that sets the scale of the theory. Here J p is an external source, Qa a 
regulating function and /a is a c-number quantity. 

The action is arbitrarily divided into a "kinetic term" and an interaction part as 



S [<f>; A] = 1 J ^- V P K X ( P 2 ) + S mt [0; A] 



where Pa stands for the regulated free propagator 



K 



P A (p») = (2*)"— i^- (1) 



n 2\ _ m~\ d 

P 

with K(z) being an arbitrary (but fixed) cutoff function which vanishes faster than any polynomial when 
p 2 — > oo and satisfies K(0) = 1. 
Using the path integral identity 

we can write the condition of the independence of the generating functional on the cutoff AdZ /dA = 
as an evolution equation for the total action S[<j), A] 

a\ / 1 f a / 2\ ( $S 8S S 2 S \ f SS ■ „ ll \ 



provided 



/*=/<*M/e- 2 M|H^-,-/^(o) 




where Qa (p 2 ) = -Pa (p 2 ) Q (p 2 ) and Q (p 2 ) satisfies Q (p 2 ) = (tj/2)Q (p 2 ). (We have defined a = da/dt, 
with t = — log(A/Ao) and Aq being a fixed scale ). In Eq (^) we have taken into account the non-trivial 
evolution of the field itself: <j) p = (ry/2)<^ p and the derivatives with respect to A should be understood as 
acting only on the coefficients of the action, not on the fields. This convention will be used throughout 
this paper. 

Note that even though Eq. (p) depends on the source J p through the Boltzmann weight, the argument 
of the v.e.v's does not. So Eq. (|) is valid for arbitrary J p if and only if the arguments of the r.h.s and 
l.h.s vacuum expectation values are equal. So from now on we will drop the v.e.v brackets. 

To complete the renormalization group transformation we should add the canonical rescalings of the 
couplings. The easiest way to proceed is to re-write Eq. (||) in terms of dimensionless quantities by 
extracting the appropriate powers of A carrying the canonical scale dependencies. We get 



• v 



S(f>pS(f)- 

• V „ 2 K '(P 2 )\ , SS f , d' 5S 



where all quantities should be understood as dimensionless. We have defined K'(z) = dK(z)jdz and the 
prime in means that the derivative does not act on the momentum conservation delta functions in 
S. This is our final ERG equation for a scalar theory. It is exact, and being supplied with the initial 
condition S\t=o defines completely the renormalization group flow of the action. 
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2.2 The action 



The second step is to define the theory, i.e. to specify the space-time on which the model is defined, 
its fields and its symmetries. Then, after characterizing the model we are instructed to take an action 
which is the most general one consistent with the stated symmetries. In our example this would be a 
power-like interaction with arbitrary number of derivatives invariant under space rotations and discrete 
^-symmetry. The next step would be to use our equation to obtain the exact renormalization flows of 
our theory. 

This program is, however, as nice as unrealistic. It is obviously impossible to deal with such a large 
action. And even in the case the /3-functions were computed somehow, the infinite set of equations would 
be intractable. Consequently, one has to develop a reasonable approximation. 

In the fermionic case we will use a mixed approximation consisting in expansion in derivatives and 
truncation in the number of fields. With this purpose in mind we will consider an example of the mixed 
approximation in the bosonic case. This will also illustrate the calculation of fixed points and their critical 
exponents within the ERG approach. A short review of other approximations used in the literature and 
discussion of their effectiveness is given at the end of the next subsection. Namely, we consider a scalar 
field theory on the Euclidean space with Z 2 symmetry and we keep all local terms up to six fields and two 
derivatives. The number of the derivatives is chosen to be the minimum one which allows, in principle, 
a non-vanishing anomalous dimension. Therefore, we take the following initial interaction action: 

(a 2 + b 2 p\) <fi Pl ■ ■ ■ P4 8(Epi) 
where oi, a 2 , 0-3, b 2 and 63 are real coupling constants. 



S int = (2tt) d / ai0 p 0_ p + (2tt) 



-3d 



(2tt) 



— 5d 



(a 3 + b 3 pf) (j> Pl 



Pi,---,Pe 



2.3 Fixed points and critical exponents 

It is a simple exercise to substitute the above action into Eq. (H|) and obtain after a bit of algebra the 
complete set of /3-functions within this approximation. The system of /3-functions will depend on the 
particular shape of the regulating function K(x) through four real quantities: 

K' (p 2 ) , [3= f p 2 K' (p 2 ) , 7 = K' (0) , 5= K" (0) . (5) 

However it can be proved that, within our approximation, the real dependence on the scheme is given 
not by four but only by two independent parameters. This simplification can be made explicit after the 
following redefinition: 

ai—>-ai, a 2 —> — a 2 , a 3 — > — =- a 3 , b 2 -> — ¥ b 2 , b 3 -> -^-^b 3 . 

7 a j or 7 cry o? z 7 z 

So, in these new variables the set of /3-functions reads, 

= V - y (126a - Sal) (6) 



0*1 


= (2 4 


rj) ai - 12a 2 


- 6b 2 /x + Aa\ 


« 2 


= (4- 


d + 2t]) a 2 - 


30a3 — 10b 3 /x + \&a\a 2 


O3 


= (6- 


2d + 3?y) a 3 - 


f 24aia 3 + 1603 


62 


= (2- 


d + 2rf) b 2 + 


\Qa\a 2 + ld>a\b 2 — 2Qb 3 




= (4- 


2d + 377) 63 ^ 


192 144 


h 


- -r-a 2 6 2 + 24ai6 3 + 24aia 3 H — r-a 2 , 





where y = S/^f 2 and x = ory/ftS. 
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The value of the coupling constants at fixed points will be sensitive to x and y. The critical exponents 
in principle are independent on them. However, the exponents depend on the scheme because of the 
truncation, as one can immediately realize noting that the anomalous dimension 77 (which in the vicinity 
of a fixed point becomes one of the critical exponents) depends almost linearly on y (some non-linear 
dependence enters through a± and b 2 ). A similar pattern of scheme dependence was found in Ref. [jl4[ . 

To be concrete, let us consider briefly the numerical results that can be obtained from Eq. (J7|) in 
d = 3. To this end we find first the fixed points of the theory: the points in the space of couplings for 
which all the /3-functions vanish. That is, we solve the following system of equations: 

= (2 + n) ai - 12a 2 - 6— + 4a? 

x 

b 3 

0=1 + 277) a 2 - 30a 3 - 10— + 16aio 2 

x 

= 377 C13 + 24ai<23 + 16a 2 

= (-1 + 277) b 2 + 16aia 2 + 16ai& 2 - 206 3 

192 144 
= (-2 + 3rj) 63 + — a 2 b 2 + 24a x 6 3 + 24 ai a 3 + —a\. (8) 
o 

where 77 is defined in Eq. (||). 

By resolving the second, third and fifth equations in (||) with respect to et 2 , 03 and 63 respectively 
one can easily reduce this system to a pair of non-linear coupled algebraic equations for a\ and b 2 . These 
remaining equations have to be solved numerically. The trivial point, where all the couplings are zero, 
is always a solution and corresponds to the Gaussian fixed point. As it is known, in the model under 
consideration there is also one non-trivial fixed point, the Wilson fixed point. The critical behavior of 
the model at this point belongs to the same universality class as the Ising model in d = 3. The system 
of equations (||) has many other zeroes besides the two corresponding to the physical fixed points. The 
appearance of spurious solutions is a common feature of the polynomial approximation, and as we will 
discuss later, is very difficult to handle (see Refs. ||, ||). However, analysing successive approximations 
it is possible to identify the physical solution among the fictitious ones: the latter do not stabilize as we 
increase the order of the approximation. In our simple scalar model we can manage high orders of the 
polynomial expansion, so we can easily find the true fixed point. This is not the case for fermionic systems 
where a six-degree polynomial in fields contains more than hundred terms within the approximation with 
two derivatives. We will discuss these issues more extensively later in this article. 

Having found the true fixed point, the next step is to study the asymptotic behavior of the coupling 
constants flow in the vicinity of it, governed by the critical exponents. They can be computed by 
linearizing the RG transformations near the fixed point. That is, if <?; is a generic coupling constant, 
then its variation in the vicinity of a fixed point g* = {g\, g%, . . . g*, . . .) is approximated by Sgi = <ji = 
Rij\ g *5gj, where <7i(t) = g* + Sgi(t) and Rij is the matrix dgi/dgj. The eigenvalues of Rij\ g * can be 
identified with critical exponents. The biggest one, Ai, is the inverse of the exponent ^, which governs 
the correlation length in the critical domains, and the second one A 2 , is minus the exponent w associated 
to the slope of the <f> A operator (5- function at the critical point. 

For our system (|7|) the matrix of linear deviations reads: 

/2-77 + 8ai -12 -6/x \ 

' 16a 2 l-2r?+16ai -30 -10 /x ' 

24a 3 32a 2 -3?? + 24ai 

16(a 2 + fe 2 ) I601 -l-2?7 + 16ai -20 

V24(a 3 + 6 3 ) 2|8 a2 + ip6 2 24ai 192/5a 2 -2 -377 + 2401/ 

where 77 takes the value given in the first line of Eq. (Q) and all the variables are evaluated at the fixed 
point. 

It is a simple exercise to carry out the above program in our example. The numerical analysis of 
the solutions of @ shows that the anomalous dimension 77 is almost linear in y (for fixed x). On the 
contrary, for any fixed y the function 77(2) has a maximum at some value x*(y). Consequently we used 



this minimal sensitivity criterion |40| to fix the parameter x at x — x*(y). Unfortunately, due to the 
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monotonous dependence of the solution on y we are unable to set it by a similar prescription. A more 
careful analysis requires the study of the dependence of the critical exponents on this last parameter, 
and, perhaps, following this direction one could fix it by applying the principle of minimal sensitivity. 
We made this study for the critical exponents v and lo but with our working precision we could not find 
any noticeable oscillation. This, perhaps indicates the poorness of the approximation. 

Within the range 0.1 < y < 1 we computed the anomalous dimension and the first two critical 
exponents [y and lo). Considering the scantiness of our truncation the results for the anomalous dimension 
and the critical exponents v and lo are quite good. As we mentioned above the values of 77 as a function 
of y are monotonous ranging from 0.016 when y = 0.1 to 0.1 when y = 1. Notice that this range includes 
the "known" value of the anomalous dimension (77 = 0.035 ± 0.005). Within our accuracy, the value of v 
in this range is constant, 

^ = 0.53 ±0.015 (9) 

and lo ranges between 

0.86 < lo < 0.96. (10) 



The known values for v and lo are, respectively, v — 0.635 ± 0.005 and lo — 0.8 ± 0.05 |41 



Now let us discuss other approximations and compare with our results. As it was already mentioned 
in the Introduction, one of the simplest approximations is the local potential approximation in which the 
interaction Lagrangian just reduces to the effective potential H (see Refs. Q - for further applications). 
In this case the flow equation is a partial differential equation, and the fixed point equation is a non- 
linear ordinary differential equation. Though this equation has a family of solutions only one of them 
is finite for every value of the field and we immediately recognize it as the Wilson fixed point For this 
approximation the anomalous dimension is zero and the values of the critical exponents depend on the 
particular ERG equation used. Using a sharp cutoff ERG equation for the effective potential the authors 
of Ref. H obtained v = 0.687, w — 0.595. In Ref. jll| these critical exponents were calculated from the 
flow equation for the Legendre effective action with a power-like additive cutoff. The results there are 
v = 0.66, w = 0.628. With a Polchinski type equation within the same approximation the authors of 
Ref. @ obtained v = 0.649 and w = 0.66. 

Further simplification of the local potential approximation can be achieved by representing the effective 
potential as a polynomial in the field || , Q . In this case the differential fixed point equation transforms 
into a system of the algebraic equations, but many spurious solutions appear. The presence of such 
fictitious solutions is thus the characteristic feature of the approximations based on truncations in the 
number of fields. Analyzing successive orders of the approximation it is possible in some cases to recognize 
the true fixed point solution, at least for not very low dimensions (for example for d > 8/3). Another and 
more serious problem is the convergence of the polynomial approximation as the order of the polynomial 
increases. Indeed there are some evidence that the polynomial approximation does not converge, at least 
for the sharp cutoff approximation H , even though the low orders give reasonably good numerical results 
for the critical exponents. Thus, for d = 3 the polynomial approximation of the Wilson fixed point with 
7 terms gives v — 0.657 and w — 0.705 §. 

In general better numerical results are obtained in the framework of the derivative expansion |]l0f - 
ifbif , where the interaction Lagrangian is expanded in powers of momenta. Since we are interested in the 
low-energy properties of the model, contributions of terms with high powers of momenta are expected to 
be sub-leading. No general solid results on the convergence of this expansion have been obtained so far 
(see, however, a discussion in Ref. fl3|]). The local potential approximation is thus the leading order of 
the derivative expansion. When the next-to-leading order (i.e. terms with two derivatives) is taken into 
account the numerical results improve in general. Thus, for the ERG equation of Ref. jll| r\ = 0.054, 
v = 0.618 and w = 0.897. For the Polchinski type equation the critical exponents turn out to be scheme 
dependent, as in our example above. In |l4j ] it was shown that for a certain reasonable variation of the 
scheme parameters the ranges of variation of the critical exponents are the following: 0.019 < r\ < 0.056, 
0.616 < v < 0.637 and 0.70 < w < 0.85. 

We see that the mixed approximation with only 6 terms, studied in this article, works rather well 
comparing to other methods. Its obvious advantage is simplicity. However, it possesses the drawbacks 
of both the polynomial approximation and the derivative expansion. We expect that corrections due 
to higher orders in the action should improve the results (^) and (|io|). In principle, they also can be 
improved by choosing the best suited renormalization scheme, as it has been done in Ref. [141. 
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What should be learned from this example is the procedure which enables to calculate the critical 
exponents starting from the expressions for the (3- functions (^). The computation for fermions within 
the mixed approximation will be basically a translation of the formulas of this section. 



3 ERG equation for fermionic theories 

In this section we derive the ERG equation for a pure fermionic field theory on the Euclidean space of 
dimension d. The equation will be very similar to that for bosons. 
The action is splitted, as usual, 

S = Skin + Sint , 

where Si nt is an arbitrary function of the fields and derivatives or momenta (if we work in the momentum 
space) and Skin is a regulated version of the usual kinetic term, 



where Pa is now the matrix 

Pa (p) = (2vr) d — . 

K (z) is a cutoff function with the same properties as that for bosons. 
The starting point of the derivation is now the relation 

^- - ^_ p P A l + x-pQX 1 ) 1 + p a^ p - Q'kXp 



tr 



(p^Pa) 5(0)) - ((V> Va 1 ~ X-pQI 1 ) p A ( p K Vp " QaXp)) 



where the trace is taken over the spinor indices. This is the counterpart of Eq. (|2j) of the previous section 
and, as there, it can be used to identify the rate of change of the kinetic term. 

Imposing the independence of the generating functional Z = < 1 > on the scale A, we come, after a 
bit of algebra, to the relation 



SSint 



+ J (|x- p QaVp + |V'-pQX 1 Xp + X-pQxVp + V'-p<3I 1 Xp) + (/a) ) 
where the anomalous dimension is similarly defined by 

A>=2^" ^p=2^V 

We can write an equation for (^s) by just adding the contribution of the kinetic term. It will be 
satisfied if 

' 6S ■ SS 8 ■ SS \ 

, tT Pa ( p > Tf TT Pa W FT - 

p \otp p Sip_ p Hp H~p) 

+ 1 [w P ( Pa {p) p ^ {p) ) +> ~ ^ ( Fa " {p) Pa {p) ) wr P 

Similarly, the equations for the terms containing Xp and Xp and the term /a without the sources and 
the fields are 

Qa (p) = Pa (p) Q (p 2 ) , /a - - J dt J Q~ 2 (p 2 ) X- p P^ (p) X P , 
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with Q (p 2 ) being a scalar function which evolves according to the equation Q (p 2 ) = (rj/2)Q (p 2 )- 

Finally, after performing the pertinent rescalings, we arrive at the equation we will use in the sequel, 

S = [2(2n) d K'(p*)(^it™-^it™) (11) 
JO [fl-d + ri n2 K'(p 2 )\{ SS , SS\ f (- u & SS , u d' 5S\ 

Notice that, once the above ERG equation is derived, it would be easy to obtain a similar one for a 
model involving also other fields and interactions just by combining the manipulations here with those of 
the previous section. (The resemblance of Eq. ( pd| ) and the Polchinski type equation for scalar theories 
is pretty evident.) 

As a final comment about our equation let us comment that we present it on Euclidean space as 
it is customary in the field. For our purposes, however, there is nothing special about the Euclidean 
formulation, as finally what one obtains is just a set of relations among coupling constants. In fact, we 
have also derived the counterpart of Eq. ( |lT| ) for Minkowski space. It is not so nice because of the presence 
of an extra imaginary unit coming from the functional derivatives of the Minkowskian "Boltzmann" factor 
e lS in the second term. Nevertheless, with this equation we have computed the /3-functions for a simplified 
action (one without operators with six fields) in much the same way we will explain later for Euclidean 
space: they are finite, real and consistent with the desired symmetries, as they should be. We have not 
proceeded further, but the parallelism between them and their Euclidean counterparts strongly supports 
the common lore that both should contain the same physical information and that the choice of space 
is much a matter of taste. Nevertheless, it would probably be nice to afford a complete calculation in 
Minkowski space. 



4 The fermionic action 

We turn now to consider spin 1/2 fields. We want to apply the machinery to a concrete model, which we 
define in this section through its symmetries, then discuss sensible truncations and, finally, comment on 
the prescriptions which we have actually used to build its action in a systematic way. 

We choose a general theory with N spin-1/2 fields on the two-dimensional Euclidean space which obeys 
the discrete symmetries of parity, charge conjugation and, to obtain reflection positive Green functions, 
reflection hermiticity (see Ref . ^| for a precise definition of them) in addition to the invariance under the 
standard Euclidean transformations. Further we will impose the invariance under the chiral symmetry 
transformations of U(N)r x U{N)l. 

The next step is to choose an appropriate truncation. A reasonable idea is to use a kind of the 
derivative expansion, which was proved to be quite efficient in bosonic theories |h| - However, 
unlike in the scalar case, the approximation without derivatives (effective potential) does not work and 
the leading order should include terms both without derivatives and with one derivative. The reason is 
simple: Eq. ( |ll| ) contains, due to the sum over polarizations, the factor j) in the Kadanoff terms, while 
a similar equation for bosons does not (see, for instance, Eq. ([|)). Another significant difference is that 
in the scalar case a general potential contains an infinite number of independent operators (powers of 
fields), whereas for finite N a general function of fermionic variables with fixed number of derivatives has 
in any case a finite number of terms due to the Grassmann nature of the field. 

However, for a given but large enough N (the case which we will be mainly interested in) the number 
of different structures grows fast as the number of derivatives increases. In practice, the approximation is 
intractable already in the order with 3 derivatives unless the degree of the polynomial of the fields is also 
restricted. We consider two approximations here. The first one is the pure derivative expansion in next- 
to-leading order (with all terms without derivatives and with one derivative) . The second approximation 
is the mixed one where the number of derivatives and also the power of fields are truncated. In this 
calculation we consider terms with up to three derivatives and up to six fermionic fields. This allows for 
non-vanishing anomalous dimension. 
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The final preparatory step is to write down the action. To construct it systematically we consider each 
symmetry in turn and derive the restrictions it imposes. For the second approximation, we will neither 
comment the details of the construction, nor write down the whole action, which in this case consists of 



107 independent operators. The interested reader is referred to |32 for a thorough discussion. Here as an 
example we present only terms S^- S ' q ^ (s stands for the number of fermionic fields and q for the number 
of derivatives) of the action with 2 and 4 fields. In the momentum representation they read: 

S (4 < 0) = gi(S 12 S 34 -P 12 P 3 4) + g2V^Vi 4 , (12) 
5 (4,2) = {mip+ 2 + m-xV\i ' P34 + m 3Pi2 2 } x (£12 £34 - P12-P34) 

+ -!>lPl2 2 + r 2Pl2 ■ P34 + r 3Pi2 2 } X ^2^34 + { s &ti Pl2° + s 2Pli ' PlZ + S3P12 P34 

+ W'pr 2 fc } x + tpiip^e* (S 12 S U ~ A2P34), (13) 

where we introduced the following notations: 

5 12 e^ fc), Pi2^T(Pi)i s r(P2), vi 2 =T(pi)i j r(p2) (14) 

(the flavour indices are summed up) and = (pk ±Pi) 3 ■ is of course the kinetic term. 

In the first case, however, the action is written as an infinite series in increasing powers of derivatives 
with coefficients being arbitrary functions of scalar operators without derivatives built out of the fermions. 
We work in the N — > 00 limit. 

From the discussion above we see that there are two independent scalar operators without derivatives, 
which in the coordinate representation are (compare with Eq. (p~4|)) 

R(x) =f (i)7 J 'f(ig''(a:)7 3 /(a;) (15) 

and b 

U(x)=ip {x)ip a (x)if) (x)ip b (x) — ip (x)j s ip a {x)ip (x)"f ' s 4> 1 '(x) 

Then, the most general action in our approximation can be written as 

S = J d 2 x{A(R(x),U(x)) 

+ ^(x)j j d j ^(x)B(R(x),U(x)) 

+ ^(x)^ j i;(x)'^(x)fd J i;(x)'^(x)f^{x)C{R(x),U(x)) (16) 

+ ^(x)j j ^(x)(T(x)d^ a (x)^ b (x)^ b (x) - {x^sdj^ix)^ {x) ls ^\x))D{R{x),U{x)) 

+ ^(x) 1 ^(x)e l ^ a {x)d l V{x)^\x) ls ^\x) - T{x)ip a (x)^ b (x) ls d.^ b (x))E(R(x), U(x)} , 

where A, B, C, D and E are arbitrary functions of the operators R and U . The standard kinetic term is 
contained in the second term of this action. Thus, the operator B should be normalized by the condition 
5(0,0) = 1. 



5 Pure derivative expansion for fermions 

In this section we will study fixed point solutions of the fermionic ERG equation in the next-to-leading 
order of the derivative expansion. The critical exponents will be analysed elsewhere. To derive the flow 
equations we write the action (^|) in the momentum representation, understanding the operator functions 
A, B, etc. as being expanded formally in powers of R and U, and then substitute it into the ERG equation 
(pr|). After long but straightforward calculation a system of equations for generating functions A(r,u), 
B(r, u), C(r,u), D(r,u) and E(r,u) can be obtained. These functions are the same as their operator 
prototypes introduced above, but with the operator arguments R(x) and U(x) being replaced by the real 
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c-numbers r and u respectively. The equations for the generating functions are the following: 



A(r, u) = 2A- 2(1 - rj)(rd r + ud u )A + AaNB - Aa(rd r + ud u )B 

+ 2a{NrC - rC - r 2 d r C - rud u C) - 2a[(r + u)D + ru(d r D + d u D)]\ 



B(r,u) 


= 7]B- 


-2{l- rj){rd r 


+ ud u )B - 


f 2 7 r(cU) 2 - 2 7U (a M A) 2 ; 


C{r,u) 


= -(2 


- 3rj)C- 2(1 


- r))(rd r -4 


-u9„)C-4 7 (9 r ^) 2 ; 


D(r,u) 


= -(2 


-3r?)D-2(l 


- 77) (r<9 r -I 


- ud u )D - ^(drA^duA); 


E(r,u) 


= -(2 


-3i7)S-2(l 


- 77) (r<9 r 4 


- ud u )E + 2 1 {d u Af + 4 7 (d r A){d a A), 



(17) 

where a and 7 are scheme parameters defined by Eq. (|^). As before the dot means the derivative with 
respect to the flow parameter t. We see that the function E(r, u) appears only in the last equation and 
is completely determined by the function A(r, u). 

There is only one solution regular in r and u at r = u = 0. It is it-independent and is given by the 
following expressions for the generating functions: 

A(r,u) = -2aN + gr, B(r, u) = 1 + 75 2 r, C(r, u) = -2 7 .g 2 , D(r,u)=0, E(r,u) = 0, (18) 

where the normalization of the kinetic term has been taken into account. The anomalous dimension 
rj = for this fixed point. Note that Eq. (jlq) is valid for any N. This solution is polynomial in r and 



is characterized by the arbitrary continuous parameter g. Since r corresponds to the operator (15) in 
the original action, the term gr in the expression for A(r, u) corresponds to the operator gV^V^, which 
gives rise to the U(l) Thirring like term. This agrees with the known fact that in the chiral Gross-Neveu 
model the sector, corresponding to this term, decouples from the rest of the model and is at the fixed 
point for any value of g. However, the coupling constant, corresponding to the SU(N) sector, is zero for 
the solution (|l8|). Because of the too low order of our approximation, we do not think that this fixed 
point can be reliably identified with any of the known fixed points of the Gross-Neveu model. 

Let us mention for completeness that the system (|T^ ) has another formal fixed point solution. The 
anomalous dimension rj — again, and the generating functions do not depend on r. The expression for 
A is given by 

. , , . _ r du' 

A(u) =A Q + Go / - — ^ , r , 

J no 1 ~ 2G ajNlii{u , /u ) 

where Aq, Gq and uq are integration constants (only two of them are independent). One can see that 
derivatives of A are not regular at u = 0. Other generating functions for this solution are expressed 
through similar integrals. It is not clear for us whether such formal solution has any physical meaning. 



6 Mixed approximation for fermions 

In this section we will consider the approximation including the terms with up to three derivatives and 
up to six field operators. Simplification due to the truncation in the number of fields allows us to increase 
the number of the derivatives in the approximation, that in its turn amounts to non-zero anomalous 
dimension, as we will explain shortly. Having constructed the corresponding action, as it was explained 
above, the next step is to substitute it into the ERG equation (|ll]) and to compute the /3-functions or 
flow equations of our model within the chosen approximation. 

One comment is in order here. In the actual calculation of the /3-functions we used an extended action 
which was larger than the one discussed in the previous section in order to have some extra check of the 
equations. Namely, we considered an action written in a basis of functionals with two fermions and one 
and three derivatives, four fermions with zero and two derivatives and six fermions with one and three 
derivatives which were restricted by U (N) symmetry, parity and Euclidean invariance only (that is, we 
imposed neither reflection hermiticity, nor charge conjugation, nor chirality). We then projected the space 
of functionals generated by this basis onto the subspace of functionals invariant under all symmetries of 
the theory and its direct complement. The required flow equations were obtained after we performed the 
first projection, while the projection onto the complementary subspace provided us a consistency check 
of the calculation. The latter defined a set of null equations, or identities, which had to be satisfied 
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along the renormalization flow. Indeed, since the initial action at t = is symmetric, appearance of any 
non-zero non-symmetric term at t > would indicate an anomaly, which, as we have argued, does not 
appear. 



6.1 Generalities 

The next step is to find the fixed point solutions, that is, the sets of values of the coupling constants 
which make all the /^-functions to vanish. They correspond to stationary points of the RG flows, thus 
providing us with the first information of how the phase diagram of the system looks like. 

S = is equivalent to the system of 106 non-linear algebraic equations. A simplification comes from 
the observation that all the coupling constants of operators with six fields must enter linearly. The reason 



is that the only source of non-linearity in Eq. (11) is the first term of its r.h.s., and neither l S^ 6,1 - ) nor <S( 6 - 3 ) 
contribute to this term within our approximation. Resolving the subsystems S^ 6 ' 1 ^ = and S^ 6,3 ' = 
with respect to the couplings of six-field operators we reduce the system to a set of only 13 non-linear 
equations, plus an equation for the anomalous dimension rj that satisfies 

t] = 4a[-TOi + m 2 + m 3 + s x + s 2 + s 3 + s 4 + 1 - 2N(r 2 + 2s 3 + s 4 )] . (19) 

Here iV is the number of flavours; a, (3, 7, S are scheme dependent parameters defined by the same 
formulas as (||); the rest of coefficients of ( |l9| ) are different coupling constants of the four-field functionals 
of our action, see Eq. ( |l3| ) . 

The appearance of the above dependence just reflects the freedom in choosing the renormalization 
scheme, or the cutoff function in our case. Furthermore, although the /3-functions contain four scheme 
parameters, we will see that, as in the scalar field case, the fixed point solution and critical exponents 
depend only on two combinations of them. 

We found the fixed point numerically for finite N, whereas for N — > 00 we could solve it analytically 
in the leading order of the 1/N expansion. Results of both calculations are presented below. 

After the fixed points are found, the behaviour of the theory near each of them is controlled by the 
critical exponents. One of them is already known once we solved our system of equations: it is the 
anomalous dimension at the fixed point. The rest are found by linearizing the ERG flow equations near a 
chosen fixed point g* in the same way as it was explained for the bosonic case in Sect. 2.3. The eigenvalues 
of Rij\ g * = {dgi/dgj) \ g * are identified with the critical exponents. They can be viewed as the anomalous 
dimensions of the corresponding operators in the vicinity of the fixed point. 

Finally, let us turn again to the issue of scheme dependence. As we have already mentioned the 
values of the coupling constants at a fixed point are scheme dependent, thus reflecting that they are 
not universal quantities. Critical exponents, on the other hand, are universal and should be scheme 
independent. Nonetheless, due to the truncation, the scheme dependence will inevitably appear. We 



will proceed in a common way 1 14 and try to find a scheme where the effect of the dependence will 
be the least. To this end, we will apply the principle of minimal sensitivity in order to fix the scheme 
dependence. 



6.2 iV^oo 

Now we are going to define the large N expansion in our model which allows to obtain analytical results 
in the leading order in 1/N. To this end we substitute each coupling constant gi by N Zi gi and study 
the limit N — > 00 keeping gi fixed, z, can be any real number, but for the sake of simplicity we restrict 
ourselves to integer values. 

We are looking for the sets {zi} for which the /^-functions for redefined couplings are finite and, if 
possible, non-trivial in the N — > 00 limit. With the above requirements, we found two different patterns 
of the 1/N expansion, which lead to inequivalent results. We will label them by I and II, and discuss in 
turn. 

The Type I solution is obtained by taking zi = —1, where i runs over all four-fermion couplings. 
With this definition, the anomalous dimension vanishes in leading order in 1/N and at the fixed point 
the coupling constant g\ of the operator (S12S34 — P12-P34) in S^ 4,0 ', Eq. (|l2|), is fixed, 

.91 = -1/(07) 
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z = 0.1 


z = 0.5 


z = 1.0 


z = 2.0 




2.258 


2.239 


2.217 


2.175 


w* 


0.122 


0.616 


1.250 


2.610 


V 


1.130 


1.128 


1.125 


1.116 



Table 1: Local minimum A| of X\(w,z), the most relevant critical exponent, for different values of z. 
w* is the value of w at which the minimum is reached and rj is the corresponding value of the anomalous 
dimension. 

whereas the coupling g2 of the operator U/2V34, giving rise to the U(l) Thirring excitations, is arbitrary. 
The characteristic polynomial -P(A), associated to the matrix Rij\ g * , can be computed analytically: 

P(A) = A 2 (A + 2) 12 (A + 4) 83 (A + 6) (A 2 + 6A - 8) 

x (-A 5 - 12A 4 + (8w - 44)A 3 + (64w - 16)A 2 + (32w + 64)A - (128™ + 256)) , 

where w = [35/ (aj) and z = 5/j 2 . The critical exponents can be read from P(X). There are 100 scheme 
independent eigenvalues, most of them coinciding with the canonical values 0, —2, —4 and —6. The non- 
trivial ones are —3 + VT7 = 1.1231... and —3 — vT7 = —7.1231.... The rest of the eigenvalues are given 
by the roots of the polynomial 

Q(A) = -A 5 - 12A 4 + (8w - 44)A 3 + (64w - 16)A 2 + (32w + 64)A - (128m; + 256) 

and are w-dependent. If w < 0, that includes, for instance, the case of the exponential cutoff function 
K(z) — e~ KZ , the most relevant critical exponent is Ai = 1.1231... 

The fixed point solution contains two free continuous parameters gi and mi. This is the expected 
result for the chiral Gross-Neveu type model because U(l) Thirring like excitations (which in our action 
are controlled by 172) decouple from the rest and this subsystem is conformal invariant (i.e. it is at the 
fixed point) for any value of 172- For the SU(N) part there exists a discrete set of fixed points, the fixed 
point of Dashen and Frishman being one of them. It is reached when the constant g\ is of the order 1/N, 
as in our case. So we can expect that our solution is the fixed point of Ref. |3q] . However, the values 
of the anomalous dimension do not match. For the cited fixed point it is non-vanishing in leading order 
in 1/N and not zero as we have found. This discrepancy with the exact result of Dashen and Frishman 
could be caused by the low accuracy of our truncation. We cannot reject, however, the possibility of 
having found a different fixed point (see Ref. |39[). 

For the Type II solution we let a combination of coupling constants to be of order 1, instead of 
0(N~ 1 ). We do not know if this possibility is an accident of the truncation or if it would survive in a 
complete calculation. The most significant feature of this solution is its non-zero anomalous dimension, 

_ 4 07 _ 4 1 w 
' 3 6a 3 6 z 

Unlike the previous case, we could not find the exact analytical expression for the characteristic 
polynomial. However, by computing numerically the eigenvalues for different values of z and w we could 
deduce some exact results. It turned out that none of the critical exponents coincide with their canonical 
counterparts. Moreover, most of them are functions of the combination ^. Thus, there are 82 eigenvalues 
A = — 8 eigenvalues equal to |(1 — ^) and 4 equal to 2 — The remaining ones are functions with 
more complicated dependence on w and z then the ratio w/z (and even a few have a non- vanishing 
imaginary part, which is not unusual in approximations based on truncations). We studied numerically 
the most relevant critical exponent, which turned out to belong to the class with non-simple dependence 
on w and z, for a wide range of values of the scheme parameters. As it also happens in the scalar case, 
for any value of z this exponent always has a minimum at some w = w* . Such behaviour suggests to use 
the minimal sensitivity criterion to fix the parameter w to its critical value w*. Unfortunately, due to 
the monotonous dependence of the exponent on the parameter z in the range analysed, we were unable 
to fix it by a similar prescription. We show in Table 1 A|(z) = Xi(w* , z) and w* for some values of z and 
the corresponding anomalous dimension. 
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6.3 Finite N 



For finite number of flavours we could not find fixed point solutions analytically, so we studied the zeroes of 
the /3-functions numerically. In general the number of different solutions of a system of coupled non-linear 
equations is not known a priori, and common programs for root-finding (such as the FindRoot routine of 
Mathematica) do not give all its solutions. Nevertheless, after some experience was acquired and relying 
on the results for high N we could choose a reasonable range of values of the coupling constants and 
examine it minutely in searching for the fixed point solutions. 

A more serious problem is to discriminate between zeros, which correspond to real fixed point solution, 
and spurious roots, which are artefacts of the truncation. This problem, which appears in the bosonic 
case too, is perhaps the Achilles' heel of the approximations based on truncations ||, ||. We present 
the class of solutions of which we are more confident. These are mainly the ones which asymptotically 
match with some solution clearly identified in the framework of the large N expansion. 

To begin with, let us select a particular scheme and find the solution for different values of N. We 
take w = —2 and z = 0.5, corresponding to the exponential regulating function K(x) — e~ x . We will 
analyse the dependence on these parameters later on. 

There is a sequence of fixed point solutions for various (interger) N which asymptotically approaches 
the type I solution of the N — > oo limit discussed above. For this sequence Nij increases with N and tends 
to 4.87..., while the most relevant critical exponent Ai decreases with TV and asymptotically approaches 
the value 1.1231..., in a full agreement with our previous results on the 1/N expansion. The second 
eigenvalue was found to be complex, that is apparently an artefact of the approximation used. Another 
mismatch is that the type of the solutions of the sequence does not fit the type of the N — > oo fixed 
point. As we wrote above, the latter is the two-parameter continuous line of solutions parametrized by 
#2, which is related to the Thirring U(l) sector, and mi, whereas the solutions for finite ./V are isolated. 
Most likely, this qualitative change of the type of the solution occurs when we go from infinite N to, 
though large, but finite N within our approximation. However, we have not studied this phenomenon in 
detail. We present in Fig. 1 the plots of Nrj and Ai as functions of N. 




N 



Figure 1: N- r\ (solid line) and Ai (dashed line) as functions of N. The corresponding fixed point solution 
matches the Type I solution of the large N limit. 

Let us discuss now the dependence of the solutions on the renormalization scheme. The parameter z 
enters the flow equations only through the anomalous dimension and, for this reason, the dependence of 
the fixed points solutions on z is rather simple: for instance, it is almost linear for ij. The dependence 
on w is more complicated, and we studied the behaviour of r\ and Ai under the change of w for fixed N 
and z. In this analysis, as in the scalar case, we looked for some non- linear w-dependence so that we 
could invoke the principle of minimal sensitivity to fix the value of this parameter and reduce the scheme 
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dependence. To this end, we take z = 0.5 and TV = 1000. The curve 77 vs. w is monotonous and decreases 
with w, while the first eigenvalue Ai reaches its minimum value Ai = 1.12511 at w = —45, which increases 
as we lower N: it is equal to 1.1273 for N = 500 (and the minimum is reached at to = —23), 1.146 for 
N = 200 (at w = -10), 1.1519 for N = 100 (at to = -8), 1.695 for N = 10 (at w = -2.4) and finally, 
2.560 for N = 3 (at w = -0.5). 

For a sequence of fixed point solutions which match the Type II solution the situation is different. 
For N large enough, (say N > 1000), the numerical solutions are in good agreement with the infinite 
TV analytical result (for example the value of r\ = 1.99 for z — 0.5 and w = —2, this is to be compared 
with the exact result r\ = 2 for A" — > oo). As we lower N the values of the anomalous dimension 77 and 
the most relevant eigenvalue Ai decrease and near N = 142 (actually at N = 142.8 if we let N to take 
non-integer values) the solution ceases to exist and joins another branch of fixed point solutions. For fixed 
points belonging to the second branch 77 and Ai remain finite as N — > 00. However, we found, that some 
couplings did not behave as powers of A" in this limit and, therefore, this solution cannot be associated 
with a fixed point in the N — > 00 limit in the sense stated previously. At the joining point rj = 1.88 and 
Ai = 5.80. We show in Fig. 2 the curves rj(N) and Xi(N). Further studies should be carried out in order 
to understand the structure of these solutions better. 



100 200 300 400 500 600 700 800 

N 

Figure 2: r\ (solid line) and Ai (dashed line) as functions of N for z = 0.5 and w = —2. In both curves 
the upper branch corresponds to the solution which matches the Type II solution of the large A^ limit. 

There are also other solutions, which were found for low N. For some of them either 77 or Ai have 
minima in w, but in other cases both curves are monotonous in w. They also display various behaviours 
as A^ — > 00. 

Finally we will consider the special case AT = 1. Fierz transformations relate U(N )-covariant local 
operators (like ip (pi) ip h (P2)) to U(N) scalars (like ip (pi) ip a (p2))- So, for A^ = 1 the Fierz transfor- 
mations establish additional relations between the available operators, which give relations between the 
coupling constants that allow us to reduce the system of the flow equations considerably. 

After a bit of algebra and discarding the trivial solution we finally get an equation for 77, 

= -120w 2 z + 288wz 2 + r](13w 2 -132wz + 210w 2 z + 288z 2 - 720wz 2 ) 
+ r) 2 (99wz - 90w 2 z - 432z 2 + 594ttjz 2 ) + t; 3 (162z 2 - 162wz 2 ). 

As in the previous analysis we had to choose some particular scheme, i.e. fix w and z, and solve the 
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= -0.1 


w = -0.5 


w = -1.0 


w = -2.0 


z = 0.1 


1.763 


3.691 


6.316 


11.747 


z = 0.5 


1.418 


1.790 


1.418 


3.388 


z = 1.0 


1.376 


1.559 


1.811 


2.349 


z = 2.0 


1.354 


1.445 


1.569 


1.834 



Table 2: Values of rj for TV = 1 for different values of the scheme parameters z and w. 

equation numerically. Unfortunately, the system of fixed point equations is not underdetermined and 
there is no room for a free parameter in the fixed point solution as it is the case in the Thirring model. In 
fact such property takes place in the previous order approximation (terms with less than three derivatives), 
however r\ vanishes identically there. This is in a complete agreement with the first solution of the pure 
derivative expansion considered in Sect. 5. The reason why this property is lost in the approximation 
with three derivatives is unclear for us. Probably, one has to go to the next order of the approximation 
to understand this issue. 

We solved the equation numerically for different values of w and z. As in previous examples the 
fixed point solutions are almost linear in the parameter z, their behaviour with w is more complicated. 
However, in the range of values studied, we did not find any non-monotonous dependence on w either of 
the critical couplings or of the anomalous dimension. We present in Table 2 some of the results for 77. 

7 Discussion and conclusions 

In this article we gave an account of main ideas and techniques of the ERG approach. We presented steps 
of the derivation of the ERG equation, which plays the crucial role, for scalar and fermionic theories. We 
discussed main approximations, which have been developed so far, for solving it and showed how some 
characteristics of models, like fixed points and critical exponents, can be computed. We also reviewed 
briefly some non-perturbative results, obtained within this approach for the scalar Z 2 -symmetric model 
in three dimensions and for the two-dimensional Gross-Neveu type model. In the latter case the results 
are quite recent. As we argue above, within the 1/N expansion the Type I fixed point solution is an 
excellent candidate for the Dashen-Frishman fixed point, whereas the other one presents evidences to be 
a new fixed point, with quite intricate properties, not discussed previously in the literature. 

We saw that in scalar theories truncation in the power of fields (like in the polynomial approximation 
of the local effective potential or in the mixed approximation) leads to numerous spurious solutions 
with the subsequent problem of identification of the true physical one, especially when the order of the 
polynomial is high enough. However, if the pure derivative expansion is used, the spurious solutions and 
the problem of identification do not appear. We demonstrated that in the case of fermionic theories the 
situation is similar. One should only keep in mind that for finite N, due to the Grassmannian nature of 
the fields, terms with I derivatives do not contain more than k max (l) = 4N+2[l/2] spinors (here [x] stands 
for the integer part of x). Thus, each term of this expansion is a polynomial in fermionic fields. When 
we use the mixed approximation, i.e. we restrict the action to the number of fields less than k max (l) for 
a given /, as it was the case in most of the examples considered in Sect. 6, we observe the appearance 
of many spurious solutions. When we carry out the pure derivative expansion, either in a general form 
(see Sect. 5) or by taking into account all possible operators with powers till k max (l) in fields, like in the 
example with N = 1 in Sect. 6.3 (I < 3 there), we see that no spurious solutions appear. Thus, for N = 1 
there are only three fixed point solutions, besides the trivial one, two of them having complex coupling 
constants, thus being rejected at once. It would be interesting to perform a pure derivative expansion 
for, say, N = 2 and check if the above features hold. 

All approximations which include derivatives contain scheme parameters, which cannot be eliminated 
by redefinition of the coupling constants. Such scheme dependence does not appear, however, in the 
lowest approximations without derivatives. This is clearly illustrated by our scalar example in Sect. 2.3. 
For the approximation without derivatives, i.e. when 62 and 63 are zero, 77 = and the equation (0) is not 
present, no scheme parameters appear in the first three equations of the system (ffl). When the terms with 
two derivatives are taken into consideration, the fixed point solution and critical exponents depend on 
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the scheme parameters x and y. This dependence on the regulating (cutoff) function is analogous to the 
renormalization scheme dependence in the perturbation theory calculations. For observable quantities, 
like critical exponents, this dependence is unphysical and is an artefact of the approximation. We used 
the principal of minimal sensitivity to "minimize" their dependence on the scheme. 

An interesting feature is the seemingly good results for the large N limit in the fermionic model. At 
the computational level this is related to the fact that at the leading order of the 1/N expansion the 
system of the flow equations simplifies dramatically and no room is left for spurious solutions. The same 
simplification occurs in scalar models as well |17j], K would be useful to understand deeper reasons 
for this. Apparently, improvement of results in the large N limit is a quite general feature of the ERG 
approach. 

Recent developments of the ERG formalism allow us to conclude that it can be regarded as a reliable 
non-perturbative method of studies of physical models. Calculational techniques, developed in scalar 
theories, were proved to be effective enough in calculations of the numerical characteristics of the models 
and analysis of renormalization flows. In fermionic theories much less work has been done so far. It would 
be interesting to probe similar calculational techniques in higher approximations, than discussed here, 
or in other physical models. For example, the ERG formalism should be extended to higher dimensions. 
The equation ( |TT1 ) is prepared for that. The number of spinor structures will increase in this case and, 
therefore, one will have to handle more terms in the action. However, we guess that, due to the greater 
complexity, two derivatives may be sufficient to obtain interesting results, or at least, according to the 
standard rule that quantum effects become less important when the dimension is increased, we hope that 
the results will be more transparent within the same approximation. 

We did not discuss gauge theories in this article. This is because of lacking of a satisfactory ERG for- 
malism for the moment. The main problem is to maintain the gauge invariance of the ERG equation and 
the renormalization flow, which is obviously broken when cutoff functions of the type (Q) are introduced 
in order to set the physical scale. Development of an effective formalism for gauge theories seems to be 
the most challenging problem in the ERG approach nowadays. 
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